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The current concordance model of cosmology is dominated by two mysterious ingredients: dark 
matter and dark energy. In this paper, we explore the possibility that, in fact, there exist two dark- 
energy components: the cosmological constant A, with equation-of-state parameter ica = —1, and a 
'missing matter' component X with wx = —2/3, which we introduce here to allow the Friedmann 
equation written in terms of conformal time n to be form-invariant under the reciprocity trans- 
formation a(ri) — > l/a(n) of the scale factor. Using recent cosmological observations, we constrain 
the present-day energy density of missing matter to be £lx,o = —0.11 ± 0.14. This is consistent 
with the standard ACDM model, but constraints on the energy densities of all the components 
are considerably broadened by the introduction of missing matter; significant relative probability 
exists even for Qx,o ~ 0.2, and so the presence of a missing matter component cannot be ruled out. 
Nonetheless, a Bayesian model selection analysis disfavours its introduction by about 1.5 log-units 
of evidence. Foregoing our requirement of form invariance of the Friedmann equation under the 
reciprocity transformation, we extend our analysis by allowing wx to be a free parameter. For this 
more generic 'double dark energy' model, we find wx = —1.02 ±0.20 and Qx,o = 0.08 ±0.57, which 
is again consistent with the standard ACDM model, although once more the posterior distribu- 
tions are sufficiently broad that the existence of a second dark-energy component cannot be ruled 
out. Moreover, the two-dimensional posterior in the (wx, f2x,o)-plane is strongly bimodal with both 
peaks offset from the standard ACDM model corresponding to (—1,0), although the latter is still 
admissible; this bimodality is in contrast to the correctly-centred unimodal posterior obtained when 
analysing simulated observations from a ACDM model. The model including the second dark energy 
component also has a similar Bayesian evidence to ACDM to within the numerical uncertainty. 



1. INTRODUCTION 

Over the past decade, cosmological observations have 
confirmed that the background expansion of the universe 
is accelerating [HH]. This remarkable phenomenon is usu- 
ally explained by assuming the existence of a single dark- 
energy component, often modelled as a perfect fluid with 
a (generally time-dependent) equation-of-state parame- 
ter w(z) that results in it exhibiting a negative pres- 
sure. The simplest form of dark energy is a cosmological 
constant A, which corresponds to a constant equation of 
state W\ = — 1. Together with cold dark matter, which is 
key to explaining the formation of structure in the uni- 
verse, the cosmological constant gives rise to the standard 
ACDM model, which provides a good fit to existing cos- 
mological observations. Nonetheless, there have been a 
large number of other exotic forms of matter proposed to 
provide alternative explanations for the current acceler- 
ating universal expansion [3j S] , including, for example, 
topological defects [5]. 

In this paper, we remain focussed on the ACDM model, 
but with the inclusion of a second dark energy compo- 
nent, which is introduced (in the first instance) to allow 
the Friedmann equation written in terms of conformal 
time t) to be form invariant under the reciprocity trans- 
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formation a(r\) — > l/a(rf) of the universal scale factor [BJ. 
Such an invariance is of general interest, but may be par- 
ticularly relevant for Penrose's recent 'Cycles of Time' 
cosmological model |7J, which posits a cyclic universe in 
which the ultimate infinitely expanded state of one phase 
(or 'aeon') is identified with the initial singularity of the 
next. 

For a homogeneous and isotropic universe described 
by the Friedmann-Robertson- Walker (FRW) metric with 
curvature parameter k, the Friedmann equation describ- 
ing the dynamical evolution of the scale factor a(t) can 
be written as 

(f) 2 = E^oa- 3 ^, (1) 

where H = a/a is the Hubble parameter (the dot de- 
notes differentiation with respect to cosmic time t), and 
the energy density pi of each of the constituent com- 
ponents of the universe is taken into account through 
the corresponding present-day density parameter fi^o = 
87rG j Oi.o/(3-ff 2 ) and its equation-of-state parameter Wi, 
which we will assume throughout to be time-independent. 

In the ACDM model, the total density parameter is 
usually taken to comprise of contributions from radia- 
tion (w — |), matter (typically modelled as dust with 
w = 0), curvature (w = —3), and the cosmological con- 
stant (w = —1). These are listed in Table |lj in which one 
can see an obvious 'gap' that we term 'missing matter' 
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curvature 




-2/3 


missing matter ? 


Ox 


-1 


cosmological constant 





TABLE I: Equation-of-state parameters for different con- 
stituents of the universe. 



with w = — |. Interestingly, forms of matter do exist for 
which w = — |, such as domain walls [8HTU] , or particular 
scalar field models [11] . More importantly, however, the 
inclusion of a component of precisely this form is required 
if the Friedmann equation written in terms of conformal 
time rj is to be form invariant under the reciprocity trans- 
formation a(i]) — > l/a(r]) [B]. This is easily seen by first 
making the change of variable drj — dt/a in the Fried- 
mann equation and including an additional missing 
matter component X; this reads 



component to the energy content of the universe, it 
is natural to extend one's investigation by allowing its 
equation-of-state parameter to vary, rather than fixing it 
to w = — | . This more generic 'double dark energy' model 
comes at the cost of breaking the reciprocity invariance 
of the Friedmann equation, but it is also of interest in its 
own right since the observed acceleration of the univer- 
sal expansion may be driven by more than just a single 
dark-energy component. 

The structure of this paper is as follows. In Section [2] 
we give a brief summary of the phenomenology of an 
additional missing matter component with w — — | by 
investigating its effect on the expansion history of the 
universe, in particular the distance-redshift relation, and 
on the evolution of perturbations, through the cosmic mi- 
crowave background (CMB) and matter power spectra. 
In Section [3j we describe our Bayesian parameter esti- 
mation and model selection analysis methodology and 
the cosmological data sets used to set constraints on our 
'missing matter' and 'double dark energy' models. The 
results of these analyses are given in Section [4] and our 
conclusions are presented in Section [5j 



1 (da 



(2) 

Making the change of variable 12(77) = l/ a (v)i one then 
immediately finds that 



PHENOMENOLOGY 
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Hi \dr, 



(3) 

Given the observed values of the various present-day 
density parameters f2j,o, this equation is clearly not 
identical to ([2]). Nonetheless, we do find that, provided 
Qx,o 7^ 0, the Friedmann equation is form invariant, 
with a swapping of roles between radiation and the 
cosmological constant, and between matter and our 
additional 'missing matter' component (hence our name 
for it). We note that the curvature density preserves its 
dynamical role under the reciprocity transformation and 
also that the remaining Einstein equations are satisfied 
for the new 'u>' values implied by ([3]), showing that the 
reciprocity transformation is a symmetry of the entire 
set of cosmological equations 1 . 



It is of interest to investigate the phenomenology of 
a cosmological model containing a second component X 
with negative pressure (in the event the energy density 
is positive), in addition to a cosmological constant. Since 
our missing matter model (for which wx = — §) is just a 
special case (albeit an important one) of our double dark 
energy model, we will focus here on the former as being 
a representative example of the latter. 

The effect of the additional component X on the global 
expansion history of the universe depends only on the 
equation-of-state parameter wx , whereas its effect on the 
evolution of perturbations will also depend on the nature 
of the component X, in particular its assumed dynamical 
properties. We therefore consider these two issues sepa- 
rately. 



Once one admits the possibility of adding an extra 



1 It is further shown in [6] that, provided Q,x 7^ 0, the Friedmann 
equation (J2| is, in fact, form invariant under the more general 
Mobius transformation a(rf) = [a — 7 a(r))]/[fi — 8 a(rj)], where o, 
0, 7 and 8 are constants, of which the reciprocity transformation 
a(?7) = 1/(1(77) i s a special case (note that there are only three 
effective constant degrees of freedom in the Mobius transforma- 
tion; a convenient choice is to set 5a — f3^f = 1). Under the general 
Mobius transformation, however, one finds that the Friedmann 



equation |2]l becomes 



1 tda 
Hi 1*7 



4 

1=0 



where Q i}0 = Qifi(a, 0, 7, 8, f2 r ,o, f2 m ,o, £lk,o> &x,o, ^A,o) are a 
new set of 'density parameters', each of which, in general, de- 
pends on all the original density parameters, as well as the pa- 
rameters in the Mobius transformation. Thus one loses the simple 
swapping of roles between different original energy densities that 
occurs for the reciprocity transformation. 
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magnitude M, namely 




m(z) = M + 51og lc 



FIG. 1: Dimensionless luminosity distance HocIl(z) as a func- 
tion of redshift z for a concordance ACDM cosmology with 
an additional component X with equation-of-state parameter 
wx = — |, for different values of Six o- 



2.1. Background expansion history 

The global expansion history of the cosmological model 
is most conveniently represented through the distance- 
redshift relation. Indeed, comparing the predicted rela- 
tion between the luminosity distance and redshift z 
of an object with observations of astronomical 'standard 
candles', such as Type la supernovae, has provided the 
most direct and convincing evidence that the expansion 
of the universe is accelerating. 

The luminosity distance to an object at redshift z is 
given by 



d L (z) = (l + z) 



y/\Qk,o\Ho 



(4) 



where Sk{x) = smhx, x, sin a; for spatial curvature pa- 
rameter k — —1, 0, +1 respectively, and the comoving 
radial coordinate x{z) is determined by the expansion 
history: 



X{z) 



dz 

W)' 



(5) 



where H (z) is obtained from the Friedmann equation ([!]). 
The inclusion of the £lx,o into Q thus directly affects the 
expansion history embodied in H(z), and hence can serve 
either to increase of decrease the luminosity distance 
4l(z) to an object at redshift z. Fig. [T] illustrates this ef- 
fect for a few representative values of flx,o- If ^x,o > 0, 
the apparent luminosity is reduced and hence the lumi- 
nosity distance is increased compared to the standard 
ACDM model. The opposite effect occurs for flx,o < 0. 

The power of the luminosity distance as a cosmological 
probe resides in the fact that it can be simply related to 
apparent brightness m(z) obtained directly from a set 
of standard candles, each (assumed to be) of absolute 



<Il{z) 
Mpc 



25, 



(6) 



where the constant offset ensures the usual convention 
that m = M for an object at dz, = 10 pc. Type IA 
supernovae constitute a set of 'standardizable candles' 
that can be used to constrain cosmological models in this 
way |12j . 



2.2. Evolution of perturbations 

An additional component X will impact on the growth 
of perturbations through its contribution to H(z) and 
the evolution of the matter density. Moreover, we assume 
here that X has the same dynamical behaviour as that 
usually assumed for a generic dark energy component. 
In particular, we use the CAMB [TJ] dark-energy mod- 
ule developed by |14j . in which dark energy is assumed 
itself to exhibit Gaussian adiabatic perturbations. It is 
worth noting that, as the equation-of-state parameter ap- 
proaches — 1, the effects of the dark energy perturbations 
disappear, as one would expect for a pure cosmological 
constant. 2 We modified the CAMB software to include 
our additional component and calculate the predicted 
power spectra of cosmic microwave background (CMB) 
anisotropies and matter perturbations, respectively, for 
several values of ilx.o! the values of the remaining cosmo- 
logical parameters were set to their standard concordance 
values. 

We plot the CMB power spectra in Fig. [2j from which 
we see that, as one might expect, the main effect of a non- 
zero f2x,o is to shift the positions of the acoustic peaks, 
which are sensitive to the spatial geometry of the uni- 
verse, and hence depend on the total energy density of 
all the components. Thus, one would expect constraints 
on ilx,o from CMB observations to be tightly correlated 
with the constraints on J1a and fife. For positive values 
of Hx,0: we also see a considerable enhancement of power 
on the largest scales from the late-time ISW effect. The 
CMB power spectrum is now well-constrained by obser- 
vations over a wide range of scales. For example, WMAP 
observations [TS] provide good constraints on the late- 
time ISW effect arising at the largest scales and on the 
first three acoustic peaks, whilst ACT data :16J accu- 
rately measures the power at high multipoles, including 
the higher acoustic peaks and damping tail, where there 
is also considerably overlap of scales with measurements 
of the matter power spectrum. 



It should be borne in mind, however, that a possible physical in- 
stantiation of an additional component X with wx = — § could 
be in the form of domain-wall topological defects, for example, 
in which case the effect on the generation and evolution of per- 
tubations may be very different to that assumed here. 
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FIG. 2: CMB power spectra for a concordance ACDM model 
with an additional component X , with equation-of-state pa- 
rameter wx = — §, for several values of flx,o- 
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FIG. 3: Matter power spectra for a concordance ACDM model 
with an additional component X , with equation-of-state pa- 
rameter wx = — §, for several values of flx,o- 

In Fig. [3j we plot the predicted matter power spectra 
for different values of fix ,oj again the other parameters 
are set to their concordance values. We see that the domi- 
nant effect of the additional component is on the normal- 
isation of the matter power spectrum. The amplitude of 
fluctuations is enhanced for flx,o > and suppressed for 
^x,o < 0. By contrast, the positions of the acoustic oscil- 
lations, which depend largely on the matter density, are 
unaffected by the introduction of the additional compo- 
nent. 



3. ANALYSIS 

We now perform a Bayesian parameter estimation and 
model comparison analysis of our 'missing matter' and 
'double dark energy' models, using recent cosmological 



Parameter 


Prior range 




[0.01,0.03] 




[0.05,0.20] 


flk,o 


[-0.05,0.05] 


e 


[1,1.1] 


T 


[0.01,0.03] 


n s 


[0.5, 1.2] 


lo gio A 


[2.5,4] 


A sz 


[0,3] 


A c 


[0,20] 


A P 


[0,30] 



TABLE II: Ranges of the uniform priors assumed on the stan- 
dard ACDM parameters in the Bayesian analysis. 



observations. In particular, we use temperature and po- 
larisation measurements from the 7-year data release of 
the Wilkinson Microwave Anisotropy Probe (WMAP; 
|15j). together with the 148 GHz measurements from the 
Atacama Cosmology Telescope (ACT; [H>]). In addition 
to CMB data, we include distance measurements of 557 
Supernovae la from the Supernova Cosmology Project 
Union 2 compilation (SCP; [T^]) and Baryon Acoustic 
Oscillation (BAO; [17 ) measurements of distance. We 
also incorporate constraints on the baryon density infor- 
mation from Big Bang Nucleosyntesis (BBN; [18]): and 
impose a Gaussian prior on the current Hubble parame- 
ter Hq, using measurements from the Hubble Space Tele- 
scope key project (HST; [19]). 

Throughout the analysis we consider purely Gaussian 
adiabatic scalar perturbations and neglect tensor con- 
tributions. We assume a modified ACDM model spec- 
ified by the following parameters: the physical baryon 
density fit, ft 2 and CDM density f^DM^ 2 , where h is the 
dimensionless Hubble parameter such that Hq = lOO/i 
kms Mpc ; the curvature density flkfi of the universe; 
6, which is 100 x the ratio of the sound horizon to angular 
diameter distance at last scattering surface; the optical 
depth t at reionisation; and the amplitude A s and spec- 
tral index n s of the primordial perturbation spectrum 
measured at the pivot scale ko = 0.002 Mpc -1 . We also 
include the Sunyaev-Zel'doich (SZ) amplitude Asz used 
by WMAP analyses, and two additional secondary pa- 
rameters incorporated in the 148-GHz ACT likelihood: 
the total Poisson power A p at I = 3000 and the ampli- 
tude of the clustered power A c . The ranges of the uniform 
priors assumed on these standard ACDM parameters are 
listed in Table [Hj Our hypothetical additional compo- 
nent is characterised by its density parameters flx.o an d 
equation-of-state parameter wx- We assume a uniform 
prior on J7a,o in the range [—1,2] throughout. For the 
missing energy model, we have wx — — §, and for the 
double dark energy model we assume the uniform prior 
wx = [-§,-§]■ 

To carry out the exploration of the parameter space, we 
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first incorporate the extra component into the standard 
cosmological equations, by performing the minor mod- 
ifications to the CAMB code [Tj5] described in Section 



2.2 



(which implement a parameterised post-Fricdmann 
(PPF) presciption for the dark energy perturbations 
[P4]). We then include into the CosmoMC software [20] 
a fully-parallelized version of the nested sampling algo- 
rithm MultiNest [5T]|H], wm ch significantly increases 
the efficiency of calculating the Bayesian evidence and 
also reliably produces posterior samples even from dis- 
tributions with multiple modes and/or pronounced de- 
generacies between parameters. A suitable guideline for 
making qualitative conclusions has been laid out by Jef- 
freys: if Bij < 1 model i should not be favoured over 
model j, 1 < B^ < 2.5 constitutes significant evidence, 
2.5 < Bij < 5 is strong evidence, while Bij > 5 would be 
considered decisive. 



RESULTS 




-0.4-0.2 0.2 




-0.4-0.2 0.2 

«X,0 



0.020.01 0.01 



For comparison purposes, we first assume no additional 
component X, in order to determine the constraints im- 
posed by the current combined data sets on the stan- 
dard ACDM model. In particular, we find the data in- 
dicate the dominance of dark energy in the form of a 
cosmological constant with Oa,o = 0.725 ± 0.015, fol- 
lowed by matter density (dark matter ± baryons) f2 m) o = 
0.275 ± 0.015, and an almost negligible spatial curvature 
Q kfi = -0.0015 ± 0.0052. We also obtain the present 
Hubble parameter Hq = 69.62 ± 1.72. The constraints 
on the other parameters {6,t, A s ,n s } remain essentially 
unaffected by the introduction below of our additional 
component X, and so we do not consider them further. 



FIG. 4: ID and 2D marginalised posterior distributions for 
density parameters in the missing matter model (note that 
f^m.o = 1— ^A.o — fifc.o — f^x.o)- The 2D constraints are plotted 
with la and 2a confidence contours. 

factor) between the missing matter model and the stan- 
dard ACDM model is B a+x ,a = -1.45 ±0.30. According 
to Jeffreys guideline [23l [24] , the inclusion of the missing 
matter component is therefore significantly disfavoured 
by current cosmological data. 



4.2. Double dark energy model 



4.1. Missing matter model 

The inclusion of a missing matter component X with 
wx = — § considerably broadens the parameter con- 
straints. In particular, we find: S1a,o — 0.837 ± 0.149, 
which constitutes an order-of-magnitudc increase in the 
error bars as compared with the standard ACDM model, 
n m ,o = 0.276 ±0.015, fJ fe>0 = -0.0047 ±0.0068 and H = 
69.90 ± 1.76. Figure [4] shows ID and 2D marginalised 
posterior distributions for the density parameters (note 
that £l m ,o — 1 — ^a,o — &k,o — ^x,o)- As expected, we 
observe a clear degeneracy between Qx,o and both Q A ,o 
and rifc.o- The ID constraint on the density parameter of 
missing matter is £lx,o = —0.11 ± 0.14. The current data 
thus prefer a slightly negative value, which is difficult 
to interpret physically, although the ID marginal shows 
moderate relative probability even for £lx,o ~ 0.2, and 
so the presence of an appreciable missing matter compo- 
nent cannot be ruled out. Our results are, however, still 
consistent with a standard ACDM model. 

This view is supported by our Bayesian model compar- 
ison. We find that the log-evidence difference (or Bayes 



We now allow for the equation-of-state parameter wx 
for our additional component to be a free parameter (al- 
beit still independent of redshift), for which we assume 
a uniform prior in the range wx = [— §,— §]■ We thus 
allow for the possibility that this second dark-energy 
component could be a form of phantom energy with 
wx < —1 [2S] ■ Figure [5] shows the resulting ID and 
2D marginalised posterior distributions for wx and the 
density parameters in the model (once again, note that 
^m,o = 1 — ^a,o — &k,o — &x,o)- At the top-right of the 
Figure we also give a representation of the 3D posterior 
in the (wx, ^x,o, ^A,o) subspace, where the colour indi- 
cates the value of il A , . 

The ID constraints on the standard parameters are 
as follows: fl A . = 0.647 ± 0.573, fl mfi = 0.275 ± 0.015, 
n k<0 = -0.0031 ± 0.0058, H = 69.95 ± 1.80. The con- 
straints on the parameters describing the additional, sec- 
ond dark-energy component may be given as wx = 
-1.02 ± 0.20 and fl x ,o = 0.080 ± 0.574, although these 
numbers obscure the bi-modal nature of the marginal 
("Wx, ^x,o)-space and (wx, f2A,o)-space distributions (see 
below). These results are clearly consistent with a stan- 
dard ACDM model, although the inclusion of the addi- 
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FIG. 5: ID and 2D marginalised posterior distributions for density parameters in the double dark energy model (note that 
^m,o = 1 — ^a.o — fifc,o — Qx.o)- The 2D constraints are plotted with la and 2a confidence contours. The top-right panel shows 
the 3D posterior distribution in the (wx, &x,o, &A,o) subspace, where the colour code indicates the value of £Ia,o- 



tional dark-energy component has again resulted in the 
uncertainties in the constraints on the standard param- 
eters being much larger than those obtained assuming a 
ACDM model. Indeed, the ID marginal for flx,o shows 
moderate relative probability even for O^.o ~ 1- 

Moreover, the 2D and 3D marginal distributions in 
Fig. [5] have some curious and interesting features that 
are worth noting. First, as might be expected, we again 
see a pronounced degeneracy between ^a,o an d Qx,o- Of 
more interest, however, is the bi- modal nature of the 2D 
marginals in the (wx, ^A,o) and («x,!lx,o) planes, both 
of which have a distinctive 'butterfly' shape, albeit ex- 
hibiting opposite correlations. Focussing on the latter, 
we see that the two peaks are offset from the standard 
ACDM model (—1,0), although that model remains ad- 
missible. Indeed, with the help of the 3D marginal in 
the (wx, &<x,o, ^A,o) subspace plotted in the top- right 
of Figure [5j we see that the two modes of the distribu- 
tion correspond to models with £Ia,o ~ 0.4, Qx.o ~ 0.3, 
wx ~ —1.2 and Qa,o ~ 1-0, ^x,o ~ —0.3, wx ~ —0.8, 
respectively. The former model has the advantage that 
the density parameters for both dark energy components 
are positive, but requires the second component X to be 
a form of phantom energy. In the latter model, the second 
dark energy component has a more physically reasonable 
value of wx , but is required to have a negative density pa- 
rameter, which is difficult to interpret physically (at least 
more so than a negative pressure). Indeed, we see that 
this latter case is broadly consistent with our findings 
for the missing matter model discussed in the previous 



subsection. 

The marginal distribution in (Qx,o, ^A,o) subspace 
shows a strong correlation between these energy densi- 
ties that would imply the potential for a trade-off be- 
tween them. One might be concerned, however, that the 
marginal distribution plotted is strongly dominated by 
the contribution (after marginalising over wx) from near 
wx = —1. If so, one could then not infer the potential of 
a trade-off between these two energy densities at (any) 
other values of wx ■ To investigate this possibility, we also 
calculated the conditional distributions in (£Ix,o>&,a,q) 
subspace for a small set of fixed wx -values in the range 
[—0.7,-1.3]. The resulting distributions were, however, 
very similar to that plotted in Figure [5j and so indicat- 
ing that the two energy densities can indeed be traded-off 
against one another. 

Also of interest is our Bayesian model comparison, 
which finds that the log-evidence difference (Bayes fac- 
tor) between the double dark energy model and standard 
ACDM is Ba+x,a = -0.16±0.30. This shows that neither 
model is preferred over the other; indeed they are indis- 
tinguishable to within the statistical uncertainty on the 
computed evidence values. Thus, the two additional pa- 
rameters flx,o and wx in the double dark energy model 
allow it the freedom to fit the data sufficiently better than 
ACDM to compensate for the corresponding increase in 
the prior volume, and hence the model is not penalised 
by the evidence. 
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5. DISCUSSION AND CONCLUSIONS 

We have investigated the possibility that there exist 
two dark-energy components in the universe: a cosmo- 
logical constant, with w = — 1; and an additional compo- 
nent X with equation-of-state parameter wx ■ In the first 
instance, we fix the equation-of-state parameter of X to 
the value wx — — §• This 'missing matter' model corre- 
sponds to the special case in which the additional com- 
ponent is required for the Friedmann equation written in 
terms of conformal time r\ to be form invariant under the 
reciprocity transformation a(r\) — > 1/0(77). Foregoing this 
requirement, we then consider the more general 'double 
dark energy' model, in which wx is a free parameter as- 
sumed to have uniform prior in the range wx — [— § , — 5] • 
For both models, we perform a Bayesian parameter esti- 
mation and model selection analysis, relative to standard 
ACDM, using recent cosmological observations of cosmic 
microwave background anisotropics, Type-IA supernovae 
and large scale-structure, together with constraints on 
the baryon density from Big Bang Nucleosynthesis and 
on the Hubble parameter from the Hubble Space Tele- 
scope key project. 

For the missing matter model, the introduction of the 
additional component X significantly broadens the con- 
straints on the standard parameters in the ACDM model, 
but leaves their best-fit values largely unchanged. The ID 
marginalised constraint on the missing matter density pa- 
rameter is Slx,o = —0.11 ±0.14. Thus, current cosmolog- 
ical observations prefer a slightly negative value, which 
is difficult to interpret physically, but the posterior on 
this parameter is sufficiently broad that significant rel- 
ative probability exits even for Qx.o ~ 0.2, and so the 
presence of a missing matter component cannot be ruled 
out. Nonetheless, our results are consistent with ACDM 
and our Bayesian model selection analysis disfavours the 
missing matter model, as compared to ACDM, by about 
1.5 log- units of evidence. 

For the double dark energy model, the constraints 
on standard ACDM parameters are again considerably 
broadened. The ID marginalised constraints on the sec- 
ond dark energy component are £lx,o = 0.080 ± 0.574 
and wx = —1.02 ± 0.20, respectively, which are again 
consistent with ACDM. Once more, however, the ID 
marginalised posterior on flx.o is sufficiently broad that 
even flx,o ~ 1-0 is not ruled out. More interestingly, 
the 2D marginal distributions in the (u>x,^A,o) an d 
(wx,&x,o) planes are both bi-modal, exhibiting a 'but- 
terfly' shape. In particular, the peaks in the (wx,^x,o)~ 



plane are offset from the ACDM value (—1,0), although 
the latter is still acceptable. The two modes of the distri- 
bution correspond to models with fi^o ~ 0.3, £lx,o ~ 
0.4, w x « -1.15 and n A ,o ~ 1-05, Q x ,o » -6.35, 
wx ~ —0.88, respectively. We also find that the dou- 
ble dark energy model has a similar Bayesian evidence to 
ACDM to within the numerical uncertainty, and hence 
neither model is preferred over the other. To investi- 
gate the significance of the observed bimodality of the 
posterior, we performed simulated observations of the 
CMB power spectrum and Type IA supernovae appar- 
ent brightness versus redshift, assuming a standard con- 
cordance ACDM model and observational data quality 
commensurate with the real CMB and supernovae data 
used in our analysis. These were combined with the same 
constraints on the baryon density from BBN |18| and 
Gaussian prior on Hq from the HST key project [19] that 
were used in the analysis of the real data. The resulting 
posterior distribution of the density parameters is illus- 
trated in Fig. [6j As one might expect, the 2D marginal 
distributions in the (wx, ^a.o) an d {wx, &x,o) planes do 
exhibit a characteristic 'plus' shape, which is indicative of 
the parameter degeneracies along the coordinate axes. In 
contrast to the bimodal posterior obtained from the real 
data, however, we see that the analysis of simulated ob- 
servations yields a unimodal posterior centred correctly 
on the input concordance parameters values. This sug- 
gests that the bimodal nature of the posterior in Fig. [5] 
may be driven by features present in the real data, but 
absent from the ACDM simulation. 

Since the double dark energy model remains viable 
with current cosmological observations, it is of interest 
to generalise it further. The bimodality we observe in 
Figure [5] might be evidence of contributions from two 
species X and Y, with wx ~ —1.2 and wy ~ —0.8. It 
might therefore be of interest to remove the requirement 
of having one fluid with w = — 1 and investigate if the 
value of either wx or wy settles near —1. Also of interest 
would be to allow the equation-of-state parameter wx (z) 
of the second dark energy component (and possibly also 
of Y) to depend on redshift We plan to investigate 
these possibilities in a future work. 
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